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p-regime of Wilson chiral perturbation theory. We discuss the range of validity of our formula 
and a possible extension to our computation in order to better understand the behaviour of the 
spectral density in a finite volume close to the threshold. 
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1. Introduction 



A theoretical description for the spectral density of the Wilson-Dirac operator is important 
mainly for two reasons. It has been recently shown in [2] that spectral observables computed with 
Wilson fermions, such as the spectral density of the Wilson operator, can be a powerful tool to de- 
termine interesting quantities as the chiral condensate. Additionally, a theoretical understanding of 
the behaviour of the spectral density close to the threshold might help in estimating stability bounds 
for HMC-like algorithms [3]. Cutoff effects and finite size effects (FSE) modify the functional form 
of the spectral density of the Dirac operator, and Wilson chiral perturbation theory (W^PT) is the 
right tool to analyze these aspects. 

In continuum chiral perturbation theory QfPT) in order to compute the spectral density pz)(/, m) 
of the Dirac operator D m =D + m with eigenvalues % + m one adds a valence quark with mass m v ; 
the discontinuity of the valence scalar quark condensate along the imaginary axis of the m v plane 
is proportional to the spectral density [4, 5]. The valence quark condensate in partially quenched 
chiral perturbation theory (PQ^PT) can be computed using the graded group method of [6] or the 
replica method of [7]. From the form of the spectral density in the continuum [8,5,2] for Nf = 2 



[pD(Y,m)] N LO = §<( 1 + 



7r|y|+m(3L 6 -l) + 

/ 2£v \ 1 ) 

(1.1) 
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one can already observe that the spectral density is a very good candidate to compute the low energy 
constant E, because the leading order (LO) expression of pp,{y,m) is directly related to £ and the 
next-to-leading order (NLO) corrections vanish in the chiral limit. 



2. Spectral density of the Hermitean Wilson operator 

For the Wilson operator Dw it is convenient to study the Hermitean operator Q m = 75 (Dw + m) = 
Q] n with real eigenvalues Giusti and Luscher [2] have shown that the spectral density pg of Q m 
renormalises multiplicatively as 

\Pq\r (^> m «) = Zpp Q (ZpX,m) , m R = Z m (m Q -m CT ) , (2.1) 

where Zp is the renormalisation constant of the pseudoscalar density. Moreover, once the action is 
improved the remaining cutoff effects are of 0(am) and they can be removed by suitable improve- 
ment coefficients. It is also useful to determine £ to define the mode number 

Mo 

N(A u A 2 ,m)=V / dXp Q (X,m), A 2 >A l >m (2.2) 

which is a renormalisation group invariant (RGI) quantity [2]. To compute the spectral density pg 
one can introduce a doublet of valence twisted mass fermions % v with untwisted mass m v = m and 
twisted mass pL v ; the spectral density is now related to the discontinuity along the imaginary axis of 
the twisted mass plane of the valence pseudoscalar condensate [9]. To compare the results obtained 
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in W^PT for pg with the continuum formula for po we recall that the two spectral densities in the 
continuum are connected by the relation pg(A,»i) = 2 po 



For all the details about our computation we refer to ref. [1]. Here we simply recall that for the 
valence and sea quark masses we have chosen to work in the /^-regime, and for the lattice spacing 
we have chosen to be in the generic small masses (GSM) regime where both masses are counted as 
0(a), viz. 

m,H v ,a~0(p 2 ), l/L,\/T~0{p). (2.3) 

The final result [1] for the spectral density of the Hermitean Wilson operator in infinite volume for 
Nf = 2 in terms if the PCAC quark mass is given by 



[p e (A,mpcAc) + Pe(-A,m P cAc)Wo = 2 \Pq(&, m PCAc)] NLO , cont 

TLX 



+ 



in 



2 

PCAC 



mi rAr A 16a 
pcac + _w 6 



(2.4) 



with a = 2Woa. There are two types of corrections to the continuum formula: an overall shift of 
0(a) parametrized by the W^PT LEC We and a modification of the shape of the spectral density 
with contributions of 0(a) and 0(a 2 ) parametrized byA = — y§- + + jji^j ■ It is important 
to notice that if we perform a non-perturbative improvement of the lattice theory the and Ws 
terms vanish, and if we improve the axial current in the determination of the PCAC mass, the Wiq 
term vanishes. We are thus left with remaining 0(a 2 ) cutoff effects proportional to Wg that go to 
zero in the chiral limit. 



3. Comparison with numerical data and range of validity 

To test our formula for the mode number (cfr. eqs. 2.2 and 2.4) we compared it with the 
numerical data published in [2]. The simulations of ref. [2] were performed with Nf = 2 non- 
perturbative 0(a) improved Wilson fermions at lattice spacing a ~ 0.08 fm and physical volume 
L ~ 2.5 fm with T = 2L. We have fixed F = 90 MeV and the renormalisation scale /I = m % \ 
we have performed a global fit at all the 3 masses available and all the values of Ar with 3 fit 
parameters: £, A and L&} From the global fit we obtain 

Z 1 / 3 = 266(7) MeV, A = -0.62(80) , L 6 = 6(l). (3.1) 

The numerical data and our global fit are shown in the left plot of fig. 1. For Z we obtain a 
perfectly consistent result with Giusti and Liischer [2] without performing any chiral extrapolation. 
We observe that we are also sensitive for a determination of while for A we obtain a value 
consistent with zero. We have also performed a fit fixing F = 80 MeV obtaining compatible results 
within errors. 

In our our paper [1] we have also compared our result with the continuum ^PT formula and 
with numerical data on a different range of A. Even if the analysis presented in ref. [1] is rather 
qualitative and more numerical data are needed, our conclusion is that W^PT describes the numer- 
ical data better than continuum #PT sufficiently away from the threshold, even if we can't exclude 

'in this proceedings with Z we denote its value renormalised in the MS scheme at a scale of 2 GeV. 
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A„[MeV] PCM J, (MeV) 



Figure 1: Left plot: result of the global fit of the data published in [2] with our formula for the renormalised 
mode number (eqs. 2.2 and 2.4). The fit parameters are E, A and Lg and we obtain % ' 2 /dof = 0.91. Right 
plot: the spectral density [pg(A,m) + pg(— X,m)]NW m the infinite volume. We used the parameters E = 
(275 MeV) 3 , ot PC ac = 26.5 MeV, F = 90 MeV, L 6 = 5, p = 139.6 MeV. The solid black line corresponds to 
the continuum %PT prediction, while the red (blue) lines correspond to the lattice W^PT prediction (0(a)- 
improved) on eq. 2.4 with a 2 W g ' = ±5 ■ 10 6 MeV 4 , corresponding to A = =f0.53. 

that adding NNLO terms in the continuum formula could lead to a comparable improvement. In 
the right plot of fig. 1 we show the NLO prediction for the spectral density in continuum ^PT (cfr. 
eq. 1.1, black curve) and W^PT (cfr. eq. 2.4) for two opposite values of A (red curve A < and 
blue curve A > 0). This plot is useful to make few observations. First we notice that in the region 
of X we have used to perform the global fit of fig. 1, i.e. with 80 MeV < X < 120 MeV, and for a 
reasonable choice of parameters, cutoff effects are very small. This is consistent with the fact that 
in the global fit we are not sensitive to determine A, that turns out to be consistent with zero. Cutoff 
effects become more visible in the range 40 MeV < X < 60 MeV and this might explain why all 
our fits in this range (cfr. figs. 3 and 4 of ref. [1]) prefer a negative value for A. In fact for A < 
we observe that over the whole range of X > 40 MeV the spectral density in W^PT is remarkably 
flat that is perfectly consistent with the striking linear behaviour of the mode number over the same 
range (cfr. fig. 3 of ref. [2]). The last remark concerns the behaviour of the spectral density close to 
the threshold: here the NLO corrections coming from the 0(a 2 ) terms become very large. On the 
other hand, we don't expect our formula eq. 2.4 to be valid close to the threshold. The reason is that 
close to the threshold (i.e. X — > m) at fixed lattice spacing a and at fixed quark mass m the relevant 
scale VX 2 — m 2 , related to the valence polar mass, might become very small. It is plausible then 
that for this scale 

1 . the GSM power counting \JX 2 — m 2 ~ ah?Q CD breaks down, and it is more appropriate to 
adopt the so called Aoki power counting (see ref. [9]); 

2. the p-regime power counting, M^L S> 1 (with Mj_ = TL\JX 2 — m 2 /F 2 ) breaks down and this 
mass scale enters the so-called £ -regime. The FSE diverge when X — > m [1]. A possible 
explanation to this behaviour is indeed that the zero-modes of the Goldstone bosons in the 
valence sector have to be treated in a non-perturbative manner as if the valence polar mass 
would be in the £-regime (for example in refs. [10, 11] both valence and sea quarks are 
considered in the £-regime). 
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4. The spectral density close to the threshold 

A possible way to overcome the second limitation illustrated at the end of the previous section 
is to adopt the so-called mixed Chiral Effective Theory [12], where some masses obey the /^-regime 
counting, and others are in the £-regime. 

We first start with some considerations in the continuum. We take N s sea quarks with mass m s in 
the /^-regime and N v = 1 valence quarks in the £-regime, viz. 

m s ~ (p 2 ), m v ~0(p 4 ), 1/L,1/T ~0(p). (4.1) 

As an intermediate step of our calculation we introduce a 0-term as follows 



F 2 r „ , + n I 



J§f 2 = - j-Tr [d^Ud^U 1 ] - -Tr u\u{x? Jl + Jt^U{x)U Q 



(4.2) 



It is important to remark that the introduction of a 0-term it is just a computational tool, and 
this remark becomes more important when we extend the framework to W^PT. For simplicity we 
consider the replica formulation 2 , where 

M = ^ = diag(wj,OT s , . . . ,m s ,m v ,m v , . . . ,m v ) , Ug = diag(eW ls 5 t r ) . (4.3) 

V * ' V v ' 

N s N T 

While in the continuum it is not important how the 0-term is introduced in the parametrization 
of the [/-field [12, 13], we have decided to introduce the 0-term only in the sea sector because it 
becomes relevant when extending our calculation to finite lattice spacing. 

It turns out [12] that among the Goldstone bosons there is one degree of freedom that becomes 
massless when we send the number of replicas to zero N r — > 0. It is the Goldstone boson that is 
a singlet under the SU{N r ) subgroup of SU(N S + N r ). To overcome this problem one treats the 
zero-modes related to this degree of freedom, labelled by t], in a non-perturbative manner [12]. 
The U -field can then be parametrized as follows 



e'lH/o U^. (4.4) 

Uo / 

In the valence sector additionally to the zero-modes field Uo, needed because the valence quarks 
are in the £-regime, we have the T]-field that makes Uo € U {N r ). In the sea sector we just have 
the additional tj -field and as a consequence the fluctuations <^ do not contain the zero-modes of the 
SU(N r ) generators nor of the t] -field. We can now shift the 0-term in the parametrization of the 
U -field as follow 

U{x)U e ^U(x). (4.5) 

The periodicity in 8 of the chiral Lagrangian allows us to write the partition function in standard 
fashion 

&(d)= V e- iv6 %, % = — / dd e m 2?(Q). (4.6) 
v^oc 27T Jo 



The same conclusions are valid if we use the supersymmetric formulation. 
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By performing an exact integration over the constant field 8 = 6 — r\ one obtains [12, 13] 

f v ^/ dU (detUo) v e^ Tr K + ^] , Zs = ms VL, (4.7) 

JU{Nr) 

from which one observes that the distribution of V is Gaussian and it is controlled by the sea quarks 
which are in the /^-regime. The computation of the spectral density is now straightforward and we 
obtain the NLO result 



P(0= E Pv(Ch£, Pv(C) = ^C[4(C)-^v+i(CVv-i(C)] J C = y£v, (4.8) 

V = — oo 

where E e g is a function of the sea quark mass m s and of the geometry of the 4-d box. 3 The 
sum over v can be easily done because we know the weight factor =^ (cfr. eq. 4.7). We observe 
that the shape of the spectral density is entirely described 'effectively' by the spectral density in 
the quenched theory, while the effects of the sea quarks are twofold. They change the absolute 
normalization by introducing an m s dependence in L e g- and they control the distribution of v, which, 
we remark, stays Gaussian only because the sea quarks are in the ^-regime. 

To extend this result to include the lattice spacing effects we need to discuss how to deal with 
the (9 -term and the 0(a) already at the LO of the chiral Lagrangian 

p2 £ r i aF 2 

JS? 2 = ^Tr [dpUdptf] - -Tr ulu(x)^ + ^U(x)U e - — Tr [U + U*] , (4.9) 
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where the mass matrix is of the same form as in eq. 4.3, but with a twisted mass (m v + ifi v z 3 ) in 
the valence sector. It becomes evident now the importance of introducing a 0-term only in the sea 
sector. Exactly as we have done for the continuum theory we redefine the £/-field as in eq. 4.5 and 
the resulting chiral Lagrangian becomes 



p2 £ aF 2 r 

JS? 2 = ^Tr [d^Ud^] - -Tr [^U f + Ujfi\ - — Tr t/^t/ + U U 



(4.10) 



The important question is to understand how to reabsorb the 0(a) term in a redefinition of the mass 
term in presence of a 0-term. One possible answer is based on a combination of introducing the 
6 term only in the sea sector and an appropriate choice of power counting. Once we have fixed 
the power counting for the masses as in eq. 4.1, where we now have mp = ■s/n^^~p^, ~ 0(p 4 ), we 
have to decide which power counting we want to adopt for the lattice spacing a. We can consider 
the following cases: 

• a ~ mp ~ p 4 (GSM-valence): at NLO, the sea quarks are effectively as in the continuum; 

• a ~ p 3 (GSM* -valence [15, 16]): at NLO, the sea quarks are effectively in the continuum; 

• a ~ m s ~ p 2 (GSM-sea/Aoki-valence): at NLO, lattice spacing corrections affect both the sea 
and the valence sector. 

It is easy to see that in the first two power countings the 0(a) can be reabsorbed in the untwisted 
valence mass m v , because in the valence sector there is no 0-term, due to the specific choice of 



See for instance [14] for an explicit expression of E e ff. 
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Uq. The 0(a) corrections in the sea sector are, on the other side, of higher order, leaving all the 
unwanted mixing between 0-term and 0(a) as NNLO corrections. 

Another important results of this mechanism is that in the GSM and GSM* -valence power count- 
ings the distribution of v is still Gaussian and the lattice spacing corrections appear at NNLO. This 
comes hardly as as a surprise because in the mentioned power countings the sea quarks, that control 
the distribution of v, are effectively in the continuum. While it is still possible that discretization 
effects will suffer from the first limitation exposed at the end of section 3, we believe that this 
framework could be the appropriate one to describe numerical data obtained with quark masses in 
the jo-regime (as usually is done with Wilson fermions). 

5. Conclusions 

We have computed the spectral density in the ^-regime of W^PT at NLO [1]. Our final result 
eq. 2.4 gives a good description of the numerical data available [2]. Lattice artefacts at NLO are 
small and vanish when m — > if we stay sufficiently distant from the threshold. Our formula can 
help both in estimating the safe region in the eigenvalue spectrum to use for the extraction of £ and 
also to give an estimate of W%. 

The theoretical understanding of the behaviour of the spectral density near the threshold is not 
complete and to fill this gap we have started a computation in the e/p-regime in PQW^PT. We 
have proposed a computational tool to deal with the simultaneous presence of a 0-term and 0(a) 
corrections in the chiral Lagrangian. 

We acknowledge useful discussions with L. Giusti and P. Hernandez. 
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